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Abstract Population synthesis is used to study the contribu- 
tion from unresolved sources to the Galactic ridge emission 
measured by EGRET. Synthesized source counts are com- 
pared with the 3rd EGRET catalogue at low and high lat- 
itudes. For pulsar-like populations, 5-10% of the emission 
> 100 MeV comes from sources below the EGRET thresh- 
old. A steeper luminosity function can increase this to 20% 
without violating EGRET source statistics. Less luminous 
populations can produce much higher values without being 
detected. Since the unresolved source spectrum is different 
from the interstellar spectrum, it could provide an explana- 
tion of the observed MeV and GeV excesses above the pre- 
dictions, and we give an explicit example of how this could 
work. 
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1 Introduction 

i 

The Galactic plane is known to be an intense emitter of 
X-rays and y-rays from keV to at least 100 GeV. At X- 
ray en ergies, it has recently been claimed dRevnivtsev et al.l 
2006) that the 2-10 keV emission can be explained entirely 
by a population of weak sources, mainly CVs. Above 50 
keV, unresolved sources also appear to be required to ex- 
plain the hard-power law emission obse rved by INTEGRAL 
dStrong et alj|2005l : iBouchet et ai]|2005h : AXPs and/or pul- 
sars are potential candidates. This raises the question of the 
situation for y-rays , in particular the range observed by 
EGRET, 30 MeV - 100 GeV. y-ray telescopes are relatively 
insensitive and reveal only the 'tip of the iceberg' of the 
sources, so most of them will go undetected unless it hap- 
pens that only strong sources exist in the Galaxy. Conven- 
tional wisdom is that the source contribution to the unre- 
solved ridge emission is at the few percent level, The prob- 
lem of the GeV excess in the diffuse emission compared to 
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the ex pected spectrum from interstellar processes dStrong et al.l 
2004a), is however a hint that the emission may have other 
components whose contribution is also energy-dependent. 
The failure to explain the 1-30 MeV emission measured by 
COMPTEL from interstellar components is another pointer 
to a source population. 

This topic has wide ranging implications: for example 
it has been claimed that the GeV exce ss is a signature of 
dark matter decay ( de Boer et all 1200511^1 , but this is only 
plausible if all viable alternatives are excluded, and source 
populations provide one important such candidate. 



2 y-ray source properties 

The 3rd EGRET Catalogue contains 271 sources, includ- 
ing 66 well-identified extragalactic sources (AGN), leaving 
about 200 for the present study. A subset of these may also 
be extragalactic. Almost the only definite Galactic identifica- 
tions are 6 pulsars (+ some other pulsar candidates) and the 
Crab nebula, which have isotropic luminosities L r (> 100 MeV) 



ranging from 10 ph s ( 10 erg s ) (Geminga) to 10 ph 
s _1 (10 35 erg s _1 ) (Crab pulsar). A summ a ry of y -ray pulsar 
observations is given by I Thom pson et alJ dl999h . The Crab 
at 2 kpc is perhaps the most luminous and distant Galac- 
tic source detected, while Geminga at 160 pc the least lu- 
minous and nearest. The Vela pulsar is intermediate in dis- 
tance (300 pc) and luminosity 10 38 ph s -1 (10 34 erg s _1 ) 
although it is the brightest source on the sky (8 10~ 6 cm~ 2 
s _1 compared to 2 10~ 6 cm -2 s -1 for the Crab). Other de- 
tected or candidate pulsars also lie in this luminosity range. 
Apart from pulsars, plausible identifications exist e.g. for 
the microquasars LS5039 and LSI +61 303, but these do 
not significantly help the present investigation. Most sources 
with Ly(>100 MeV)< 10 33 erg s _1 are invisible to us with 
present instrumentation. The data are too sparse to construct 
a luminosity function directly, so we leave this as a parame- 
terized input. 

1 but see critique by Bergstrom et al. (2006). 
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We approach the problem using population synthesis, 
comparing the models both with source counts and the inten- 
sity of diffuse emission^. Comparison in various sky regions 
increases the discrimination power of the comparisons. Thus 
source counts in high-latitude regions constrain the low- 
luminosity, nearby source population while counts in low- 
latitude regions constrain the high-luminosity sources, and 
diffuse emission in low-latitudes constrains both low and 
high-luminosity populations. 

As usual in such studies a major uncertainty is the spa- 
tial distribution of sources, but with plausible assumptions 
this can be modelled. Once the spatial distribution is as- 
sumed, the observed (source-produced) gamma-ray sky de- 
pends only on the luminosity function. 

For the truly diffus e (interstellar) emiss ion we make use 
of the galprop models (Stro ng et al.ll2004al) . 

The present study is designed to be independent of the 
physical details of the sources, making use of, for exam- 
ple, theoretical y-ray luminosity functions of pulsars only 
for guidance. We do not want to restrict attention to known 
classes of objects. The essential input is only geometry, the 
inverse square law, simple luminosity functions, the EGRET 
source catalogue and EGRET skymap data. 



3 Previous studies 

Most previous population studies have been aimed at de- 
ducing the nature of the EGRET unidentified sources, either 
generally or in particular as pulsars. Although the nature of 
the sources is not our goal here, these studies are still rele- 
vant. A population synthesis study of the 2nd EGRET cat- 
alogue dKanbach et al.1 IT996) showed that the unidentified 
sources had luminosities >100 MeV in the range 6 1 34 — 
3 1 35 erg s" 1 , with 700-3400 objects in the Galaxy. Flux 
dis tributions of EGRET unident ified sources were constru cted 
bv lReimer & Thompsonl d200lh and lGehrels eTaT] d2000l) in 
attempts to deduce the properties of unidentified sources, 
including division by variability and angular distribution. 
IChen et al (2001) made a population synthesis study of the 
3EG catalogue, including disk and isotropic source distribu- 
tions. 

Other studies considered pulsars specifically. Very de 



and also studied the properties of 38 low-latitude unidenti- 
fied sources, proposing an a correlation wit h SNR and OB 
associations. Yadigarogl u & Romanil d 19971) proposed that 
the EGRET unidentified s ources are compatib l e with young 
pulsars, while in contrast iBhattacharya et a l. (2003) claim 
they are better traced by spiral arms and molecular clouds. 

In view of the uncertainty in the nature of the unidenti- 
fied sources, a flexible approach to modelling is desirable, as 
described in this paper. 



4 Population synthesis 

A general-purpose population synthesis code has been writ- 
ten. Sources are assigned a density p(R,z,Ly) and sampled 
by standard Monte-Carlo techniques. Oversampling is used 
to reduce statistical fluctuations. The density is normalized 
to p at R — 8.5 kpc, in units of sources kpc~ 3 . Power-law 
luminosity functions within given Ly limits can be gener- 
ated. The (R,z) source distribution is here based on pulsars 
dLorimeril2004l) as representative of y-ray sources, but other 
distributions are also possible and will be addressed in fu- 
ture work. The resulting source list is analyzed to generate 
differential source counts N(S) and the total emission spec- 
trum both above and below a given detection threshold. The 
total spectrum is then combined with inte rstellar emission 
models from galprop dStrong et alj|2004al) . We use source 
counts for one energy range (> 100 MeV) only; it would 
be preferable to consider the energy-dependence of N(S) but 
the spectral information in the available catalogues is lim- 
ited. In future (e.g. for GLAST) this will be feasible. A note 
about beamed sources, in particular pulsars. For the present 
purpose a populations of randomly-oriented beamed sources 
is fully equivalent to a population of unbeamed sources with 
a lower spatial density of sources having an isotropic emis- 
sion. Therefore we do not explicitly include beaming in our 
population synthesis. 



4. 1 Pulsar luminosity function 

We use pulsars just as an guide to the choice of luminosity 



tailed p hys ical pulsar population synthesis studies by lGonthier (ftfflij tion. For this we use the luminosity as a function of spin- 
( 200l) and lGonthier et alj 42004) using their polar-cap model down power E : Ly oc E^ for which a wide spread exists in the 
were aimed at establishing the relation between radio and literature depending on the model. The luminosity function 
y-ray properties a nd predicting detections for future mis- can be then estimated as follows: N(L Y ) = = % §■ = 

sionslike GLAST.Earding & Zhan^M^ Gonthier^L_^ * dE x * ^ for constant birthrate Since £ocrVA 
(2003h: lHarding et al.l(l2004l) : lGonthier et al.ld2005hlCheng et alF IS^dEdT tor constant birthrate. Since Eoc B 

d2004l) considered the possibility that gamma-ray sources in P/P 3 , dE/dt oc E/P 2 and hence N( Ly) oc B~ 1 Ly 
the Gould's Belt are nearby pulsars. 



(l+2j3)/2j3 



According to the polar-cap model of Gonthie r et al. ( 2002 ) 

IZhang & Cheng! d 19981) proposed that pulsar s can account j3 « j, so N(L r ) °*L~ 2 The slot-gap model o flMuslimov & Harding] 



for the diffuse GeV excess. |Zhang et al.| 42OO0D made a de- J2003t) gives ]8 = \ so N{L y ) °c L.; 3 .The outer-gap model of 
tailed population synthesis based on the outer-gap model, ^mgjXj^ ^ gives j8 = 0.4 - 1 so N{L y ) « L" 2 - 3 - 

Ly 15 . The dynamic range of y-ray luminosity in these mod- 
els is about 1000. 



2 We will often refer to unresolved emission generically as 'diffuse' 
regardless of its true nature 
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4.2 Known populations 

We start with the EGRET catalogue sources (excluding AGN 
identifications), attempt to reproduce their source counts by 
population synthesis and hence estimate the contribution of 
unseen members of this population to the diffuse emission. 
The luminosity function is assumed to be a power law, the 
index and limits being free parameters. In our reference model 
the range of luminosities considered is Ly(>100 MeV) = 
covering the range of detected pulsars as 
discussed in the Introduction. The local density p is fixed 
by the requirement that the low-latitude source counts are 
reproduced. 

For EGRET we use a limiting flux (> 100 MeV) of 10" 7 
cm -2 s _1 ; fainter sources are detected by EGRET at high 
latitudes but in the plane the limit is higher: 1.6 10~ 7 \b\ < 
10°, 0.7 10~ 7 \b\ > 10°). For a luminosity function index - 
1.5, p = 37 kpc -3 (model lb), the low-latitude source counts 
in 300° < / < 60°, \b\ < 10° (hereafter region H) are re- 
produced (37 above threshold), the high-latitude sources are 
very few (so they must be extragalactic in this case), and 
about 6% of the Galactic emission ( 2 10~ 4 cm -2 s -1 ) in this 
region comes from the 4000 sources below the threshold. 
The fluxes above and below the threshold are about equal, 
1.5, 1.2 10" 5 cm -2 s" 1 >100 MeV in the simulation, com- 
pared to 1.4 10~ 5 cm -2 s _1 for the 37 EGRET sources in 
region H. 

For a luminosity function index -2.0, model lc, again 
choosing p so that the low-latitude source counts in region 
H are reproduced, 13% of the Galactic emission comes from 
the 26000 sources below the threshold. 

How critical is the luminosity function shape ? Flatten- 
ing the luminosity function index from -1.5 to -1.0, about 
4% of the Galactic emission comes from the 728 sources 
below the threshold. Steepening the luminosity function in- 
dex to -2.5, about 28% of the Galactic emission comes from 
the 10 5 sources below the threshold. So even a major differ- 
ence in assumed luminosity function shape does not change 
the conclusion that a significant contribution to the diffuse 
emission must come from sources physically like those in 
the EGRET catalogue but below the detection threshold. A 
steeper luminosity function or a lower minimum luminosity 
leads to a larger contribution from unresolved sources. 

The high-latitude source counts are not very constrain- 
ing, due to the low space density for these high luminos- 
ity sources, but the predicted counts are consistent with (i.e. 
do not exceed) the observed counts. For a steep luminosity 
function (index -2.0) or a low minimum luminosity (10 35 ph 
s _1 ), a significant fraction of the high-latitude sources can 
be Galactic. With the small statistics a quantitative compari- 
son of observed and predicted shape of N(S) is difficult, any 
index considered here reproduces the observed counts rea- 
sonably, and no distinction between the various indices can 
be made. 



4.3 Unseen / dim populations 

We turn now to source populations mainly below the detec- 
tion threshold, with the aim of placing constraints on their 
properties. The relevant models are 2, 3 and 4 in Table 1 . 

The luminosity function is assumed to be a power law 
with index -1.5 and various luminosity ranges. The normal- 
ization is chosen to give about a tenth of the EGRET dif- 
fuse emission, as an illustration of a plausible level. Con- 
sider first L 7 (>100 MeV) = 10 36 - 10 37 ph s" 1 , model 2. 
For p = 500, 15% of the diffuse emission in region H is 
from the 5 10 4 sources below the threshold and 44 of the 
EGRET unidentified high-latitude sources are Galactic and 
17 (cf 3rd EGRET Catalogue: 37) sources are above the 
limit in region H. Longitude profiles (not shown here) indi- 
cate that the sources below the threshold cause fluctuations 
in the emission which however are not distinguishable from 
the interstellar emission. For L r (>100 MeV) = 10 35 - 10 36 
ph s" 1 , p = 2 10 3 , (model 3), 7% of the diffuse emission 
in region H is from the 2 10 5 sources below the threshold, 2 
sources are above threshold in region H, and 19 at high lati- 
tudes. ForL 7 (>100 MeV) = 10 34 - 10 35 ph s" 1 , p = 2 10 4 , 
(model 4), 7% of the diffuse emission in region H is from the 
2 10 6 sources below the threshold, no sources are above the 
threshold in region H, and only 5 at high latitudes. It follows 
that for L 7 (>100 MeV) < 10 35 ph s' 1 all of the 'diffuse' 
emission could come from sources without violating source 
counts anywhere on the sky (e.g. by scaling up p in model 
4). This is of course highly unlikely given our knowledge of 
interstellar emission processes but cannot be excluded from 
y-ray data alone. 

If we consider that all the high-latitude unidentified sources 
are Galactic then a dim dense population like in models 2,3,4 
is necessary in addition to the bright population. Strict limits 
are then set by the requirement of not violating the observed 
diffuse emission considering plausible interstellar emission. 
If instead they are extragalactic, the dim Galactic population 
is not required. 



5 Comparison with physical pulsar population synthesis 

It is interesting to see how our simple generic approach matches 
detailed specific models. In their pulsar population synt hesis 
based on their polar cap model, iGonthier et al.l d2004l) find 
26 pulsars detectable by EGRET, which presumably would 
mean a substantial fraction of the unidentified sources are 
pulsars. They predict that 600 pulsars will be detectable by 
GLAST for a threshold 2-5 10~ 9 cm" 2 s" 1 . This matches 
best our model lb (index -1.5) 

The outer gap model of Zhang et al. d2000h predicts 32 
pulsars detectable by EGRET, which again would mean a 
substantial fraction of the unidentified sources are pulsars. 
This model predicts 1 180 GLAST pulsars, for assumed thresh- 
old of 4 10 -9 cm -2 s _1 matching best our model lc with 
luminosity function index -2.0. 
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Fig. 1 Differential source counts in 300° < / < 60°, \b\ < 10° (region H). Blue, connected points: models la-d, luminosity index -1.0,-1.5,-2,-2.5 
(left to right, top to bottom). Red unconnected points: 3rd EGRET Catalogue, excluding AGN identifications. 



Table 1 Summary of population synthesis models. Using thresholds Secret = 1 10 7 cm 2 s 1 . Sclast = 4 10 9 cm 2 s 1 . For comparison, 
the EGRET measured diffuse emission in region H (300° < I < 60°, \b\ < 10° ) is 2 10~ 4 cm~ 2 s _1 
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6 Sources can produce the MeV, GeV excesses 

The unresolved source fraction is energy-dependent, so that 
one can ask whether it can produce the well-known GeV ex- 
cess over the standard cosmic -ray interaction models. Con- 
sider first model lc (Fig 3c) with spectral index -2.0 and 
luminosity function index -2.0, combin ed with the 'conven- 
tional' interstellar emission model from Stro ng et"ail 12 004a). 
but with the cosmic -ray source distribution from lStrong et al.1 



(2004b). The GeV excess is not reproduced. However for 
this source spectrum the source contribution to the COMP- 
TEL (1-30 MeV) and INTEGRAL (.02 - 1 MeV) ranges 
might provide an explanation of the excess above the in- 
terstellar emission at those energies. Sources wit h a Crab- 
pulsa r -like index of -2. 1 w ith a break above 4 GeV dFierro et alJ 
1998; Ku iper et "all 1200 lh would be too steep to reproduce 
the GeV excess, but the COMPTEL diffuse emission could 
be fitted. 
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Fig. 2 Spectra in region H (300° < / < 60°, \b\ < 10° ) for models la-d. luminosity index -1.0,-1.5,-2,-2.5 (left to right, top to bottom). Sources 
below (dotted, cyan) and above (da shed, cyan) the EGR ET detection limit are also shown together with sources below the limit added to the 
conventional interstellar model from Strong et al. (2004a) (continuous, cyan). Data: EGRET, COMPTEL. 



galdef ID 49 5001 80 
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Fig. 3 Sources like Geminga, Vela pulsars can produce the GeV excess: spectra in 330° < I < 30°, \b\ < 5°(region A of lStrong et all C2004a) ) 
forL(>100 MeV) = 10 36 - 10 39 s _1 . luminosity index 2.0, spectral index -1.5, break at 2 GeV to -2.0 to match Geminga, Vela, p = 150. Sources 
below (dotted, cyan) and above (da shed, cyan) the EGR ET detection limit are also shown together with sources below the limit added to the 
conventional interstellar model from Strong et al. (2004a) (continuous, cyan). Data: EGRET, COMPTEL. 
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Consider now sources with hard spectra like the Geminga 
(index -1.42, 3 0-2000 Me V) and Vela (index -1.62, 30-1000 
GeV) pulsars (iFierro et al.lll998l) (also B1706, B1055 have 
hard spectra with a GeV break); we adopt a spectral index - 
1 .5, with a break at 2 GeV to -2.0, luminosity index 2.0, p = 
30 (cf. model lc); the GeV excess is easily produced (Fig 
[3J. Indeed these pulsars show maxima in E 2 F(E) around 2 
GeV very reminiscent of the Galactic GeV exce ss. This ex- ~~ 
planat ion of the GeV excess was also proposed bv lZhang & Che5 ffi erences 
d 19981) on the basis of their outer-gap pulsar model, but our 
result is not dependent on a particular physical model. 

However the diffuse emission below 30 MeV is not ex- 
plained by such sources. It is then tempting to propose a 
mixture of Crab-like and Geminga/Vela-like pulsars to pro- 
duce both the MeV and GeV excesses and reproduce the 
entire ridge spectrum. 

The spectra of the unidentified lo w-latitude sources s cat- 
ters broadly around -2 (-1.7 to -2.7) dZhang et al.ll2000h but 



6. A combination of source populations combined with the 
conventional model of interstellar emission could explain 
the full COMPTEL/EGRET Galactic ridge spectrum. 

7. Whether GLAST can settle these issues depends critically 
on the source luminosity function. 
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7 Will GLAST resolve the issue ? 

Although GLAST will not detect all the y-ray sources in the 
Galaxy, it will resolve essentially all the source flux of the 
bright populations. Thus for model lc (luminosity index - 
2.0), 75% of the source flux is above a threshold 4 10~ 9 
cm" 2 s , compared to 35% for the EGRET threshold. For 
model lb (luminosity index -1.5), 96% of the source flux 
is above the GLAST threshold, compared to 55% for the 
EGRET threshold. For the dimmer populations (L y (>100 
MeV) < 10 37 ph s _1 ) progressively less of the flux will 
be resolved. Even for L 7 (>100 MeV) = 10 36 - 10 37 ph s _1 
(model 2) only 44% of the flux is above the GLAST thresh- 
old, so the analysis will remain a challenge. 



8 Conclusions 

1 . Modelling the contribution from unresolved sources is es- 
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2. The contribution from unresolved sources to the EGRET 
low-latitude emission is at least 5-10%, and can be 20% for 
steep luminosity functions. 
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without violating EGRET source counts anywhere on the 
sky. 

4. The GeV excess can be produced naturally by a sufficient 
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5. Crab-like sources can produce the COMPTEL excess but 
not the GeV excess. 
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